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Abstract 

Using the general argument in Borel resummation of perturbation theory that hnks the divergent 
perturbation theory to the nonperturbative effect we argue that the nonperturbative effect associ- 
ated with the perturbation theory should have a branch cut only along the positive real axis in the 
complex coupling plane. The component in the weak coupling expansion of the nonperturbative 
amplitude that gives rise to the branch cut can be calculated in principle from the perturbation 
theory combined with the exactly calculable properties of the nonperturbative effect. The real- 
ization of this mechanism is demonstrated in the double well potential and the two-dimensional 
0{N) nonlinear sigma model. In these models the leading term in the weak coupling expansion of 
the nonperturbative effect can be obtained with a good accuracy from the first terms of perturba- 
tion theory. Applying this mechanism to the infrared renormalon induced nonperturbative effect 
in QCD, we suggest some of the QCD condensate effects can be calculated in principle from the 
perturbation theory. 
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I. INTRODUCTION 



The usual perturbation in the weak couphng constant in field theory is an asymptotic ex- 
pansion. When the perturbation series is sign alternating it may be resummed, for example, 
in the manner of Borel resummation. However, when the series is not sign alternating, it 
usually implies the presence of a genuine nonperturbative effect, and the Borel resummation 
of the perturbation series alone is in principle not sufficient for an adequate description of 
the true amphtude. 

The nonperturbative effect and the perturbation expansion are not totally independent; 
the former controls the large order behavior of the latter. Even with this relation, however, 
explicit calculations show that the nonperturbative effect cannot be calculated from the 
perturbation series even when the latter is known to all orders. 

The purpose of this paper is to argue, by taking a closer look at the general but heuristic 
argument that relates the nonperturbative effect with the perturbation theory, that some 
parts of the nonperturbative effect, which usually include the leading piece in weak coupling 
expansion, can be calculated in principle from the Borel resummation of perturbation theory 
along with the exactly calculable properties of the nonperturbative effect. 



II. CALCULABLE COMPONENT IN NONPERTURBATIVE EFFECT 

A general but heuristic argument that relates the nonperturbative effect with perturbation 
theory goes as follows.^ Let A{a) be an amplitude with perturbation expansion in the 
coupling constant a: 

oo 

Aia)=Y.ana-^' (1) 



and assume that A{a) be real for a > 0. We shall further assume that a„ at large orders is 
nonalternating in sign. In general, a„ diverges factorially due to renormalons or instanton- 
anti-instanton pairs, and in principle Eq. (|1|) is meaningless unless some kind of resummation 
is performed on the divergent series. To do a resummation, consider a new series 

oo 

i?(a)=5:a„(-ir+ia"+i (2) 



which is obtained from Eq. (|l|) by flipping the sign of the coupling constant. Since for 
a > this series is alternating in sign, it can now be Borel resummed, yielding a resummed 
amplitude i?PT(a), 

EpT(a) = / dhe~-B{h) (3) 
Jo 

where 

oo |'_iNn+l 

B{h) = Y. ^ (4) 


One would expect that A{a) could be obtained from the Borel resummed -BpT(a) by ana- 
lytic continuation from the positive real axis to the negative real axis in the complex a-plane. 
The problem is, however, 5pT(a) is expected to have a branch cut along the negative axis. 



^ A good introduction can be found in [0, D . 
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and consequently -Bpt(— «±ie) will have an imaginary part. Therefore analytic continuation 
alone of -BpT(a) cannot reproduce A{a) which is by definition real for a > 0. The expected 
resolution of this problem is that the true amplitude A{a) has a nonperturbative amplitude 
^Np(«) in addition to the BpT^{—a), so when they are added together the imaginary parts 
from each amplitude cancel each other, rendering the total amplitude to be real. That is, 



A{a) = ApT{a) + Aj^p{a) 



(5) 



with 



Im [Apt {a ± ie)] + Im [A^p (a ± ie)] = , (6) 

where Ap^ia) = Bp^^—a), which may be called the perturbative amplitude, has now a 
branch cut along the positive real axis in the a-plane. Performing the analytic continuation 
in Eq. we obtain 



where^ 



ApT(a ± ie) 



Mb) 



ooitie 



Oitie 



dbe~^A{b) 



(7) 



-B{-b) 
°° a 



(8) 



Note that, as we make analytic continuation of -BpT(a) from the positive real axis to the 
negative real axis in the a-plane counterclockwise (clockwise) in the upper (lower) half plane, 
the integration contour also should rotate counterclockwise (clockwise), hence the ± sign in 
Eq. (0). 

Using the dispersion relation on the resummed ApT{(y) and Eq. (^), we have 



ApT(a) 



1 /•°° Im[ApT(a' + ie)] 



da 



TT JO 
-| oo 



TT 



n=0 



a — a 

°° lm[A^T{a' + ie)] 



a 



in+2 



a 



n+l 



Thus the perturbative coefficient in Eq. (|) can be written as 

1 \m[A^p{a + ie)] 



TT JO 



a 



n+2 



da 



(9) 



(10) 



which makes the relation between the nonperturbative effect and the perturbation theory 
explicit. At large values of n the dominant contribution in (jT^) comes from the small a 
region, and so for large order behavior only the weak coupling limit of the nonperturbative 
effect is required. 

Now with Eqs. (^ and (P) the amplitude A{a) for a > can be written as the sum of 
the real parts of the perturbative and the nonperturbative terms: 



A{a) = Re [ApT(a ± ie)] + Re [ANp(a ± ie)] . 



(11) 



^ Throughout the article whenever a Borel transform is defined through a perturbation series hke (H) it 
is assumed that the value of the Borel transform at a point beyond the convergence disk is obtained by 
analytic continuation. 
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Equation (|^) also shows that the imaginary part of the nonperturbative amphtude A-^p{a) 
can be calculated in principle from the Borel resummation of the perturbation theory. The 
real part, however, is in general not calculable from the perturbation theory. 

The argument hitherto is well known, perhaps except for Eq. (J^) which allows us to relate 
the imaginary part from the analytic continuation in the a— plane to that arising from the 
Borel integral. Now, our observation, which will play a crucial role throughout the paper, is 
that Eq. (|^) suggests A^p{a) have a branch cut along the positive real axis in the complex 
coupling plane, in order to cancel the imaginary part coming from the branch cut in the 
Borel resummed perturbative amplitude. This rather straightforward observation can have 
an important consequence: it renders some part of the nonperturbative amplitude to be 
calculable from the perturbation theory.^ In the weak coupling limit [a — > 0) the approxi- 
mate functional form of A^p{a) can be rather easily determined from other nonperturbative 
techniques such as renormalization group argument or instanton calculations, and the com- 
ponent, typically the leading one in weak coupling expansion, that could give rise to a branch 
cut along the positive real axis can be easily identified. We then further fix this component 
to a more specific form by demanding that it give a branch cut only along the positive real 
axis; We expect the nonperturbative effect should not have other branch cut, for instance, 
such as one along the negative real axis, since it would imply that the perturbation series 
(0) is not Borel-summable even for a < 0. The constraint on the functional form from this 
step turns out to be sufficient enough for us to relate the real part of the above-mentioned 
component to its imaginary part, consequently rendering the former to be calculable from 
the Borel resummation of the perturbation theory through Eq. (^. 

In the next two sections we consider a few definite examples and show how this procedure 
can be realized in model calculations. In these examples we shall focus on the nonperturba- 
tive effect associated with the first singularity on the positive real axis in the Borel plane. 



III. THE DOUBLE WELL POTENTIAL 



The quantum mechanical double well potential with the action 



S 



dt (12) 



has instanton solutions, and the nonperturbative effects due to the instantons give rise to 
singularities on the positive real axis on the Borel plane, causing same sign perturbation 
series. 

Consider, for example, E{a) defined by 

E{a) = ]^[E^{a) + E^{a)-l\ (13) 

where a = \^ and i?_(a), -E+(a) are, respectively, the energies of the ground and the first 
excited states. The Borel transform 

°° a. 



E{h) = E (14) 



n=0 



^ The idea of rebuilding nonperturbative effects from perturbation theory is an old story. To the author's 
knowledge it was first speculated in 0], and an approach similar to ours was observed in H]. 
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of the perturbation series for E{a), 



E(a) = £a„a"+\ (15) 



n=0 



is expected to have multi-instanton caused singularities at 6 = 2nSo, where n = 1, 2, 3, • • •, 
and 5*0 = 1/6 is the one-instanton action. 

The large order behavior of the perturbation ( plSf ) is controlled by the first singularity at 
b = 1/3. The nonperturbative effect that causes this singularity is due to the contributions 
of instanton-anti-instanton pairs, and can be calculated from the potential of an instanton- 
anti-instanton pair, and reads 0] 



i?Np(«) = — e-V3« 



In ( -^^ + 7e) (^1 - y a) ~ y " + ^("^ 1^ 



(16) 



where 7e is the Euler constant. 

Notice that i?Np(tt) has a branch cut along the positive real axis in the a-plane, in 
agreement with our argument in the previous section. The minus sign in the logarithmic 
term, which causes the branch cut, arises from the required sign flip in the coupling constant 
to pick up the nonperturbative effect from the attractive potential of the instanton-anti- 
instanton pair 0. 

The imaginary and the real parts of E^p{a) now read 



Im[ENp(a±ie)] = i-e^^/^" 



1 + —a + O(a^) 
6 



(17) 



a 

Re[E^p{a ± ie)] = — 6"^/=^" [- ln(a) + ln(2) + 7e + 0(a)] . (18) 

Note that the real part has terms, for example, the constant terms within the bracket in 
Eq. (0), that have nothing to do with the imaginary part. These terms represent genuine 
nonperturbative effect and cannot be calculated from Borel resummation of perturbation 
theory. 

Since in this example the nonperturbative effect can be calculated in weak coupling 
expansion there is no real need to attempt to extract the nonperturbative effect from the 
perturbation theory. However, for the sake of argument, let us suppose that we knew only 
that the nonperturbative effect in weak coupling expansion was given in the form: 

E^p{a) oc — e"^/^" [ln(a) + subleading terms] . (19) 
a 

In any event, inferring this form may not be so difficult since the preexponential factor 1/a 
can be obtained by counting the number of (quasi)zero modes, in this case two, and the 
logarithmic term arises from the instanton-anti-instanton potential at large distance. 

We can now improve the form of the nonperturbative effect ([T9|) by demanding that it 
have a branch cut only along the positive real axis in the a— plane. Because the branch cut 
can arise only from the logarithmic term we can immediately see that Ej<jp{a) must assume 
the following form: 

C 

E^p{a) = e"^/^° [In(-a) + subleading terms] (20) 



5 



with C an unknown real constant. Of course, a comparison with Eq. ([T6| ) shows that the true 
value of the constant must be C = 1. We now show that the leading term, the logarithmic 
term, in the real part ( ]T8| ) can be calculated from the perturbation theory starting from the 
ansatz (pOl). The imaginary and the real parts from this expression (^0]) are then 

C 

Im[^Np(a ± ie)] = ±— e"^/^" [1 + subleading terms] (21) 

a 

C 

Ref^Npfa ± ie)] = e"^/^" [ln(a) + subleading terms] . (22) 

na 

To determine the leading term in the real part we now only need to fix the constant C. This 
constant becomes the residue, up to a calculable normalization, of the first singularity in 
the Borel plane, and can be calculated in perturbation theory |^ ^. In fact, for the Borel 
resummation 

EpT(a ± ie) = / e--E{h) dh (23) 

J 0±ie 



to have imaginary parts that can cancel the imaginary terms in Eq. (^) the Borel transform 
E{h) must have a singularity at 6 = 1/3 of the form 



9C 
7r(1^6)^ 



m = -37^^11 + 0(1 - 36)] . (24) 



To determine the residue we now consider a function R{b) defined by 

R{b) = E{b){l-3b)\ (25) 

The difference between R{b) and E{b) is that the former has a much softer singularity. 
Although it may appear that R{b) is regular at 6 = 1/3, it is easy to see that that is not the 
case. In fact, for the imaginary parts from the Borel integral ( p3D to cancel the imaginary 
parts in Eq. (|17|), E{b) should have an expansion around the singularity 
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m- 



n{l-3by 



^3 

1 - — (1 - 3b) + 0[(1 - 3bf ln(l - 36)] 

lo 



(26) 



R{b) is, therefore, logarithmically [multiplied by (1 — 36)^] singular at 6 = 1/3, but bounded. 
In terms of R{b) the constant is given by 



9 V3 

The essential point for the perturbative calculation of the residue is that the right-hand side 
of Eq. (|27|) can be written as a convergent series. The perturbation series for R{b), 

oo 

R{b) = Y.rJ)- (28) 

n=0 

is convergent on the disk |6| < 1/3 (note the boundary is included). Being bounded, though 
singular, at b = 1/3, R{b) can be evaluated in series at b = 1/3. 

We can now do some numerical checks to see how rapidly the series (pSf) ,when evaluated 
at 6 = 1/3, converges to the known exact value. From the perturbative coefficients for E{a) 
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TABLE I: Sum of the first + 1 terms of the perturbation series for the normalized residue 



oo — C'oo — 1 



N 





1 


2 


3 


4 


5 


Cn 


0.349 


0.175 


0.339 


0.487 


0.631 


0.759 


Cn 


0.349 


0.109 


0.502 


0.650 


0.862 


0.994 



given in ||^ the coefficients Vn in Eq. ( pq ) can be obtained. In Table |I| we give the first terms 
of Cn defined by 

C.^ --!: ..(-) . (29) 

Note that Coo = 1- The numbers show that from Ci to C5 the series approaches the true 
value in a steady pattern. 

Since the positions of the singularities for E{b) are known we can improve the convergence 
using an "optimal" conformal mapping 

w = w{b) . (30) 

An optimal mapping for our case is that R{b{w)) become as smooth as possible within the 
convergence disk of the perturbation expansion in the w— plane 

00 

R{biw)) = J2rnW''- (31) 

n=0 

An obvious strategy for an optimal mapping is to push all other singularities in the Borel 
plane except for the first one far away from the origin. Here we consider a mapping 



1-Jl- 36/2 
w = ^, ' . (32) 



1 + ^1 - 36/2 
This maps the first singularity to w = wq, where 

Wo = ^ ~ 0.171 (33) 

and all other singularities to the unit circle. Because the singularities (other than the first 
one) in the w— plane are relatively farther away from the origin than in the 6— plane, we 
expect the series in the w— plane to give better convergence. In fact, the first terms shown 
in Table | of C*n, which is defined by 

N 

C^ = -^Y.^n<. (34) 

^ n=0 

show a definite improvement in convergence (note again Coo = !)• 

The result of this exercise shows that the leading term in the nonperturbative effect 
caused by the instanton-anti-instanton pairs on the ground state energy can be calculated 
accurately (99% accuracy) with only the first six terms of the perturbation series for the 
ground state energy. 
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IV. THE TWO-DIMENSIONAL 0{N) NONLINEAR SIGMA MODEL 



This model is exactly solvable in expansion and mimics many interesting features 
of quantum chromodynamics (QCD). It has an asymptotic freedom, dimensional transmu- 
tation, and most interestingly for us the infrared (IR) and ultraviolet (UV) renormalons 
at the next leading order in expansion. Moreover, it was through the studies of this 

KM that the link between operator product expansion (OPE) and IR 



11, 12, 13 



model 

renormalons suggested by Parisi [1^ has become more transparent. This model can also 
provide a nontrivial test of our proposed mechanism. 

The two-dimensional 0{N) nonlinear sigma model is defined by the (Euclidean) action 



^=2 



1 f ^ 



N 



a X cr ix] — 



(35) 



where a{x) is an auxiliary field, and / is the coupling constant. At the leading order in 1/A^ 
the a fields get dynamical mass 

= ;,2e-i//(A«) ^ (36) 
where /x is the renormalization scale, through the vacuum condensate of the auxiliary field 

(0|a(0) |0) = -ViVm^ (37) 



Equation (36) also defines the renormalization group (RG) running of the coupling constant. 
The (3 function in the leading order in 1/A^ is therefore given by 



Pif) = /i' 



(38) 



We shall now test our proposed mechanism with the truncated two-point function of the 
a fields 



(39) 



which is known to all orders in OPE at order 1/A^ via the exact calculation of the self energy 
ZI(p^) in [ITH]. r(p^) can be expanded in powers of m^/p'^, which corresponds to an OPE, as 



np' 



p 



Co{m) + Cr{f{p))— + Oi — 
P \P\ 



(40) 



We keep here only the first two terms because we are focused on the nonperturbative effect 
associated with the first IR renormalon. The first term contains the usual perturbation 
expansion, and the second term, which comes from the vacuum condensate of a{x), is the 
nonperturbative amplitude that gives rise to the first IR renormalon. The terms of higher 
powers in m? /p^ are associated with the higher renormalon singularities, and shall be ignored. 
At the leading order in 1/A^, Cq = Ci = 1. At order 1/A^ they have a rich structure and 



read 16 



ooitie 



1 /-ooil 

Co(/(p)±i6) = - / 

N Jo±ie 

C^{f{p)±ie) = 



db 



-b/fiv) I 4-FW(6) + G(°)(6)l - ffW(6) 



ooitie 



N Jo± 



db 



-b/fiP) 



JiP) 



F«(6) + 6-^(6) -H'^^\b) 



(41) 
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where 



= 1 

= 1 + 5,(6), (42) 



and 



F^'){b) = b'-l 

G^^\b) = -1 + 262 -46+1 + (1-62) [^(1 + 6) +^(2 -6)+ 27s] 

BQ{b),Bi{b), whose exact forms are not important for us, are analytic functions. 

Note the renormalon poles at 6 = n, where n is a nonzero integer, in G^^\ H^^\ To avoid 
these poles the integration contour can be either on the upper or the lower half plane; for 
consistency, however, an identical contour should be taken for Cq and Ci. 

We shall identify the first term in Eq. (^Op as the perturbative amplitude Fpx and the 
second term as the nonperturbative amplitude Fnp. Performing the integration over 6 in 
Eq. ( piD it is then easy to see that Fnp at order 1/A^ is given by 

2 

T7? 

NTMfip) ± ie) = [ln/(p) + ivr + real const. + 0(/)] (44) 

= _e-i//(p) [ln(-/(p) + ie) + real const. + 0(/)] , (45) 

where in the last step we have used Eq. (^). Thus to the leading order in weak coupling 

iVrNp(/) = -e-i/^ln(-/), (46) 

which has a branch cut along the positive real axis in the coupling plane, again in agreement 
with our proposed mechanism. The origin of the imaginary part in Eq. ( ^41) lies with 
the ambiguity in obtaining the renormalized condensate, in this case the condensate of 
the auxiliary field a{x), from the dimensionally regularized condensate in 2 + e dimension 
0, [TT|, |12[ . When the perturbative and the nonperturbative amplitudes are added together 
this imaginary part is canceled by the imaginary part in Fpt(/(p) ± ie) coming from the 
pole at 6 = 1 in G(°)(6). 

From Eq. (^Tj) the perturbation expansion for Fpt(/(]5)) can be easily obtained 

oo 

rPT(/(p)) = In fip) + const. + a„/(p)"+i , (47) 

n=0 

with 

ao = -2, 

an = n\ [(1 + (-1)") ({n + 1) - 2] (for n > 1) , (48) 
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where ( denotes the Riemann C— function. The logarithmic term in Eq. (^) arises from the 
anomalous dimension of the a fields. 

It is instructive to see now that the logarithmic dependence in ( PI ) can be obtained 
without knowing the exact solution for the a self-energy. To see this let us make the 
renormalization scale dependence of F explicit, which was hitherto implicitly suppressed.^ 
r(p^,/i^, f{fi)) can be expanded in OPE as 

= p'[Coi^^Vp^fi^^)) + ^,c,i^^'/p^fi^^ma\o)^^ 

-|-higher dimension terms] , (49) 

where Ci are the Wilson coefficients, and the ignored terms involve operators of dimension 
four or higher. From the RG equations for the coefficient Ci(/i^/p^, /(/i)) and the condensate 
(0|a|0) the second term in (^9]) can be easily written up to a momentum independent factor 
as 

rfip) 

exp - / [27.(/')+7a(/')]//3(/')^/' , 



Ci(l,/(p))e-^/^(^) 



(50) 



where ■jaif) and 7a(/) denote the anomalous dimensions for the a and a fields, respectively. 
Comparing this term with the second one in Eq. (^OD we can identify 



Ci(/(p))=Ci(l,/(p))exp 



f{p) 



(51) 



Therefore the logarithmic term in Eq. ( |1^ ) can arise from a logarithmic term at order 
1/iV in the Wilson coefficient C*i(l,/(p)), and the anomalous dimensions 7{o-,q}, which are 
nonvanishing at order 1/N: 

7K«}(/)oci^[/ + 0(/^)]. (52) 

What is noteworthy here is that the contributions from these sources when added together 
conspire to absorb the imaginary parts in Eq. ( p4|) into a cut-function. If, for example, 
the logarithmic term had a different coefficient than that given in Eq. (|4^) , then after the 
absorption of the imaginary parts into a cut-function [ln(— /)] there would have remained a 
logarithmic piece (ln(/)) which has a cut along the negative real axis in the coupling plane. 
This would have failed our argument that the nonperturbative amplitude has a cut only 
along the positive real axis. 

Now, once we know the form of the anomalous dimensions for the a and a fields, and 
the logarithmic dependence in the Wilson coefficient at order the nonperturbative 

amplitude can be written with no help from the exact solution for the a self energy as 

iVrNp(/(p)) oc e-^/^(P)[ln f{p) + subleading terms] . (53) 

Then demanding rNp(/) have a branch cut along the positive real axis in the coupling plane 
we can refine Eq. (^) to 

NT^p{f{p)) = Ce-^^^^P^[\n{-f{p)) + subleading terms] , (54) 



This was allowable because the renormalization scale dependence in r(p^) can be factored out. 
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with C an undetermined real constant, from which we can obtain a relation valid in leading 
order in weak coupling: 

Re[rNp(/)] = T-Im[rNp(/±ie)]ln/(p) 

TT 

= ±-Im[rpT(/±ie)]ln/(p). (55) 

TT 

Because the imaginary part of the perturbative amplitude can be calculated from Borel 
resummation this relation renders the real part of the nonperturbative amplitude to be 
calculable from the perturbation theory. 

The leading term in rpT(/) depends only on the constant C. We can calculate this 
constant from perturbation theory using the method used in the previous section for a 
similar purpose. First, we note that for the rpT(/) to cancel the imaginary part coming 
from Eq. at positive f{p) the Borel transform of rpT(/) should have the following 

singularity at 6 = 1 (the first IR renormalon): 

fpT(6) = Y^[l + 0(l-6)] . (56) 

This shows that the constant C becomes the residue of the renormalon singularity and can 
be written as [0, |^ 

C = R{1) (57) 

with 

i?(6) = (l-6)fpT(6). (58) 

Because of the UV renormalon at b = —l the residue C cannot be directly evaluated by 
the perturbation expansion of R{b) around the origin. To map away other renormalons than 
the first IR renormalon we introduce a conformal mapping 



vm - ji - b/2 

w = ' (59) 

VTTb + J I - 6/2 



which sends the first IR renormalon to if = 1/3 and all other renormalons to the unit circle. 
Because in the if —plane the first IR renormalon is the closest singularity to the origin, C 
can now be evaluated by plugging w = 1/3 into the following series expansion of R{b{w)):^ 

oo 

R{b{w)) = Y.r^w^. (60) 

n=0 



The first terms of Cat, which are defined by 

N .1 

^^ = E^nU) ' (61) 



n=0 

were calculated using the perturbative coefficients (^8]). The numbers in Table show that 
the residue can be determined with good accuracy from the first terms of the perturbation 
theory. 



^ More precisely, it is the expansion of R{b) of the Borel transform for [rpT(/) — In / — const.] in Eq. ( |47| ) 
rather than that of the Borel transform for rpT(/). Either of the Borel transforms can be used because 
both have an identical renormalon singularity at 5 = 1 . 
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TABLE II: Sum of the first A*" + 1 terms of the perturbation series for the renormalon residue 

{Coo = -1). 



N 





1 


2 


3 


4 


5 


— Cn 


2.000 


2.000 


0.100 


0.945 


0.916 


1.052 



V. THE QCD CONDENSATE EFFECTS 

We now come to the potentially most interesting application of our proposed mecha- 
nism. Because of the nonperturbative nature of the QCD vacuum, operator condensates ap- 
pear ubiquitously in low energy QCD phenomenology, especially in the Shifman-Vainshtein- 
Zakharov (SVZ) sum rule formalism [jl^ |18[. The effects of these condensates become 
stronger at lower energies and become phenomenologically more important. They are in 
general not calculable, and treated as free parameters to be fitted by experimental data. 
The condensates are generally introduced through OPE. In Borel resummation of the QCD 
perturbation theory they appear as the nonperturbative amplitudes which are required to 
cancel the imaginary parts arising from the IR renormalon singularities. The purpose of this 
section is to see the implication of our proposed mechanism on the nonperturbative effects 
caused by these condensates. 

In general the form of a nonperturbative effect due to the condensates can be determined 
by OPE and the RG equations for the associated Wilson coefficients and the condensates. 
Once the form is determined then we can further refine it by demanding the nonperturbative 
amplitude have a branch cut only along the positive real axis in the coupling plane. This 
will then allow us to write the real part of the amplitude in terms of the imaginary part that 
can be calculated from Borel resummation. 

Of course, one should remember that not all nonperturbative effects in QCD could be 
related to perturbation theory. It is clear, for example, that perturbation theory cannot 
have any bearing on the nonperturbative effects arising in chirality violating processes. 

For definiteness, we shall consider the Adler D-function in massless QCD defined by 



D(Q^) = -W^^^^^-1, (62) 



where n(g^) (g^ = ~Q^) is the vacuum polarization function in the Euclidean region {q"^ < 0) 
of the current j'^(x) = u{x)''j^ d{x) , with u^d denoting the up and down quarks. D{Q'^) can 
be expanded in OPE as 

-(-higher dimension terms . (63) 

As before we focus on the nonperturbative effect associated with the closest singularity 
to the origin on the positive real axis in the Borel plane, in this case the first IR renor- 
malon, and ignore terms of dimension six or higher since they are associated with the higher 
IR renormalons. (OIO4IO) is the renormalization scale invariant gluon condensate of the 
dimension-four operator 



a. 



(64) 
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where G^^, is the gluon field strength tensor, and as is the strong couphng constant. [3{as) 
is the QCD /3— function: 



-a. 



(3q + (3ias + 0{al 



(65) 



where for N,. colors and Nt quark flavors 



/5o 



1 /II 



47r V 3 



34 



(4vr) 



-AT 



10 



/ 



-iVciV/ 



(66) 



The Wilson coefficients Cq, C4 can be expanded in power series in the strong coupling as{Q)'- 

00 



n=0 

27r2 



i + 5:.^f«.(g) 



(68) 



n=l 



where the coefficients dl^\w!^'^ are real numbers. 

We shall now identify the first term in Eq. ( |63D as the perturbative amphtude DpT{as{Q)) 
and the second term of gluon condensate as the nonperturbative amplitude D^p{as{Q)). 
DpT can be expressed as a Borel resummation of the perturbation series ( p7| ) 



DpTiasiQ) ± ie) = 



where Z)pt(&) is defined by 



(69) 



(70) 



This series is expected to have a finite radius of convergence (|6| = 1) set by the UV renor- 
malon at b = —1. Beyond the radius of convergence D{b) is assumed to be obtained by 
analytic continuation. Dpx is now expected to have an imaginary part with sign ambiguity 
due to the first IR renormalon at 6 = 2 in the Borel integral. This imaginary part is to be 
canceled by the imaginary part from D^p. Now -Dnp, the second term in Eq. (|53D, can be 
written as 

DMc^siQ)) a,(Q)-''e-2/'3«"»(«)C'4(«.(Q), ) (71) 

via the weak coupling expansion of 

(OIO4IO) 



oc exp 



-2 



MQ) da' 



(3{a') 



(72) 



which comes from the RG invariance of the gluon condensate. Here, 



u = 2A//5, 



' 



(73) 
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which is noninteger for most of the combinations of Nc,Nf. C4, which comes from C4 and 
the weak couphng expansion of Eq. (^), is real and calculable in perturbation, and can be 
expanded in power series:^ 

C'4K(Q),A,«^f ) = l + w^'^as{Q)+wPasiQY + --- , (74) 

where w^^ depends on and wf'\ i < n, respectively, in Eqs. (|65|) and (|68D . 

Following the argument in Sec. |I| we demand the imaginary part in Z^np come from a 
branch cut along the positive real axis in the coupling plane. Presence of UV renormalons, 
which gives rise to a sign alternating large order behavior, does not affect this requirement 
, since UV renormalons can be mapped away using a conformal mapping from the Borel 



integration contour, thus causing no essential problem for Borel resummation [15, M. Since 



in the nonperturbative amplitude ( fflD a branch cut can arise only from the factor ^^(Q)^^, 
with u a noninteger number, we conjecture that the nonperturbative amplitude is given in 
the form: 

DMc^siQ)) = C [-«.(g)]~"e-2/^°°»(«)C'4(a.(Q), A, ) , (75) 

where C is an undetermined, dimensionless, real constant. In this specific case the power 
of the coupling constant in the preexponential factor depends on u only. But, in general 
the power depends not only on u but also on the one-loop anomalous dimensions of the 
associated operators. In such a case, we propose that the branch cut arises likewise from the 
preexponential power term in the coupling, and the correct form for the nonperturbative 
amplitude can be obtained by flipping the sign of the coupling constant in the preexponential 
factor. 

The argument leading to Eq. (^) shows that when u takes an integer value for a particular 
combination of Nc and Nf, the nonperturbative amplitude Dnp cannot have an imaginary 
part. This implies disappearance of the first IR renormalon singularity in the Borel plane 
for an integer u. A further comment on this point will follow shortly. 

Now with Eq. ( [75| ) we have the real and the imaginary parts, 

Re[DNp(«(Q) ± ie)] = C cos(z/7r)a,(g)-'^e-2/^°"^(^)C'4, (76) 
Im[DNp(a((5) ± ie)] = ±C sin(z/7r)a,(Q)-'^e-2/'3o°=(«)C'4 , (77) 

from which a relation between the real and the imaginary parts is obtained: 

Re[DNp(«s(Q) ± ie)] = ± cot(z/7r)Im[DNp(a.(Q) ± ie)] ■ (78) 

This has an important implication. It relates the usually incalculable real part of the non- 
perturbative amplitude to its imaginary part that is calculable from Borel resummation. 
Moreover, this relation is not for some part only of the amplitude as in the previous two 
examples, but holds to all orders in perturbative expansion of C4. This implies that as 
far as the gluon condensate effect is concerned the nonperturbative effect can be calculated 
completely from the Borel resummation of the perturbation theory. 



^ When the higher renormalons are taken into account the series expansion for C4 is an asymptotic expan- 
sion, and should be Borel resummed. However, since we ignore all higher renormalons C4 is assumed to 
be well defined by the series expansion. The error on C4 by this assumption is 0(e~^/'^°"('3)) which is 
due to the dimension six condensates. 
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From Eq. ( ffSD and the fact that the imaginary parts in DpT(as(Q) ±ie) and Dnp(«s(Q) ± 
ie) cancel each other, we can write the Adler function in terms of Dp^ only: 

L)(q2) = [j^g ^ cot(z/7r)Im] DpT{as{Q) ± ie) . (79) 

Thus both the real and the imaginary parts of the Borel resummation are required to rebuild 
the true amplitude from the perturbation theory. The imaginary part comes from the region 
beyond the first IR renormalon {b > 2), and for its calculation analytic continuation to 
the region beyond the radius of convergence of the perturbative Borel transform ( [7DD is 
required. A more convenient method, though equivalent to the analytic continuation, is to 
use a conformal mapping to map, for example, all the renormalon singularities, or all the 
renormalon singularities except for the first IR renormalon, to the unit circle. A conformal 
mapping of the first kind was used in rebuilding the imaginary part of a met ast able- vacuum 



energy in a quantum mechanical model ||20|, and the second kind was recently used in Borel 



resummation of the real part of the Adler function which may also be used for the 
calculation of the imaginary part. Here, instead, we shall try to evaluate the constant C 
which would give a rough estimate of the nonperturbative amplitude D^p{as{Q))- 

As in the previous examples, this constant becomes the residue, up to a calculable con- 
stant, of the Borel transform at the renormalon singularity at 6 = 2. In fact, for the Borel 
resummation (^) to have an imaginary part that can cancel the imaginary part (|77D the 
Borel transform Dpx(&) should have the following singularity at 6 = 2: 

DpTib) = ^^(/5o/2)i+'^(l- 6/2)-^-^ [1 + 0(1 -6/2)] 

-|- Analytic part , (80) 

where "Analytic part" denotes terms that are analytic at 6 = 2. In general the analytic 
part cannot be calculated. Note that, as previously mentioned, the renormalon singularity 
disappears when u takes an integer value. This is obvious for a negative integer u, which 
happens, for example, at Nc = '2,Nf = 8 with u = —10, or at Nc = 3, A^/ = 15 with 
u = —176. The disappearance of the renormalon singularity for the latter case was noticed 



before |2^. What seems to have been unexpected is that the singularity also disappears 
for a non- negative integer for example, at Nc = Q,Nf = 12 with u = 1. In this case 
the singularity disappears because of the vanishing residue. The residue vanishes since the 
constant C should always be bounded, for the nonperturbative amplitude ( fTSD cannot be 
divergent. 

With Eq. (|D|), C can now be obtained by [0, § 

C = Ti-u)i2/P,y+^Ri2) (81) 

with 

R{b) = (1 - b/2y+^DpT{b) . (82) 

As in the previous examples R{2) may be evaluated as a perturbation series. However, it 
is unlikely to obtain a good estimate of the constant C by directly following the procedures in 
the previous examples, since too few perturbative coefficients are known. Only the first three 
terms of the perturbation series for the Adler function are presently known. To improve the 
situation we will exploit the renormalization scale independence of the Adler function. To 
do this, we replace in Eqs. (p^)-(pD|), and ( [7^ ) the couphng a{Q) with the running coupling 
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'^iCQ)j with ^ defined by the renormahzation scale yU^ = C^Q^- Then the perturbative 
coefficients in these equations should change accordingly as d^^ dn{0^ ''^n^ ~^ ''^niO^ 

^ WniO, with dn{^ = 1) = Wn{^ = 1) = and Wn{^ = 1) = Also the Borel 
transform ([70| ) should be redefined as 

Z^„(,0^£MlfAV. (83) 



n=0 



n! \/5o^ 



Now the renormahzation scale invariance of the gluon condensate allows us to rewrite Eq. 
(M) as 

^Np(«.(eQ)) = [-«.(eQ)]"'^e-2/«°"=(««)C'4(«.(eQ), , (84) 

and consequently the renormalon singularity corresponding to Eq. (^) is given by 

DpT{b,0 = ^(/5o/2)^+'^(l-6/2)-i-ni + 0(l-6/2)] 

+Analytic part . (85) 

Note that C in Eqs. ( PH ) and (^) is the same one defined in Eq. (^), and is independent 
of the scale ^ but dependent on the renormahzation scheme. 
With Eq. (|85|) , C can be written as 

C=ir(-z/)(2//3o)^+'^i?(6 = 2,0 (86) 

with 

R{b,0 = {l-b/2Y+''DMb,0- (87) 

We now proceed to evaluate R{b = 2,C,) as a perturbation series. Because of the UV renor- 
malon at b=—l the evaluation point 6 = 2 is beyond the convergence radius of the series for 
R{b,C,) around the origin, hence we have to map away all other renormalons except for the 
ffist IR renormalon, using a conformal mapping like 



VT+b - Jl - 6/3 

^= , V^ ^' (88) 

VTTb +Jl- 6/3 



which maps the ffist IR renormalon to w = l/2 and all other renormalons to the unit circle. 
Now in the w— plane C can be expressed as a convergent power series, that is, C = Coo where 
Cat is defined by 



1„, ../I 



with r„(^) coming from the expansion 



n=0 



W^),0 = E^n(0^"- (90) 
n=0 

Although C is ^-independent, in general Cn will have a ^-dependence because of its finite 
order summation. This property is generic for any finite order QCD perturbation series and 
can be used to improve the convergence of the series by demanding that at an optimal ^ 
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FIG. 1: Renormalon residue vs renormalization scale 



the scale dependence of the series be minimal p3[. Applying this idea to our problem we 
can hope that a better estimate of the constant C can be achieved by taking CAr(^^) at an 
optimal ^0 at which the unphysical ^-dependence disappears locally, 



0. 



(91) 



Using the calculated next-next-leading order Adler function pj, ^ and the estimated 
0(q;^) coefficient |21|, in the MS scheme at Nf = 3 quark flavors = 1)^ 



D(q2) = ^(Q) _^ X g398 ^(Q)2 _^ g 3-^xQ ^^g^3 ^ 4o)a(Q)4 + o{a^) , (92) 

where a{Q) = as{Q)/TT, we give the last two calculable terms for Cn in the MS scheme: 

C2(1.7) = 5.37, C3(2.2) = 6.96. (93) 
In this calculation we took the estimated value d^^^ = 25, which is from the recent estimate 



using a technique called "bilocal expansion of Borel amplitude" |2T[. This value is also in 
consistency with the well-known estimate in |26|. Note that the optimal ^ for Cn is at 
^ 1.7 for C2 and ~ 2.2 for C3 (see Fig. [l|). 

Because C is evaluated with the perturbation series (pO|) at w = 1/2 which is on the 
boundary of the convergence disk. Cat in Eq. (p3D should be regarded only as a rough 



^ dniO in terms of dl"'' can be found in pTI 
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estimate. The speed of convergence of the sequence Cn is expected to be sensitive on the 
size of the analytic part in Eq. (|85|), since the size of the singular term of R{b,C,) at b = 2 is 
determined by this analytic part. 

Using this estimate of C we can now evaluate the real part of the nonperturbative am- 
plitude. From Eqs. (^) 



Re[DNp(as(Q)±ie)] ^ ^3(2.2) cos(i^7r)a,(Q)-'^e-2//3o°4Q)c'. 



4 



C,{2.2)cos{urry-^ [l+w?^as{Q)+0{al)\ (94) 



where Wi^^ is defined in Eq. (|68|) and we have substituted C3(2.2) for C, Aj^ denotes the 
MS renormalization scale. 

Now this nonperturbative amplitude can be translated to a gluon condensate in the OPE 
(|63|). We define the gluon condensate by this nonperturbative amplitude by 

{^G^G^nNP = ^Re[DNp(«s(Q)±ie)] 
- ^C3(2.2)cos(z/7r)A| 



0.005 GeV^ (95) 



where we have used Aj^ ^370 MeV for Nf = 3 quark flavors . 



One should not compare the estimated value ( p5| ) directly with the phenomenologically 
fitted gluon condensate, for example, from the QCD sum rule. In the QCD sum rule, the dif- 
ference between the Borel resummed perturbative amplitude and the sum of the first terms 
in the perturbation series is approximated by power corrections, and therefore the phe- 
nomenologically fitted condensate includes contributions not only from the nonperturbative 
amplitude but also from the perturbative amplitude. One may try to extract the gluon 
condensate effect in the resummed perturbative amplitude through Eq. (|69|), for example, 
by computing the minimal term of the perturbation series with the large order behavior 
given by the renormalon singularity (|80|) . We believe, however, that this is not necessary, 
and also not a good way to handle the renormalon effect. A better approach to incorporate 
the renormalon effect, with only the first few terms of the perturbation series available, is 
to write using Eq. (|82D the Borel transform Z)pt(6) in the Borel integral (|69|) as 



:i-b/2y 



DpTib) = 7^^^ (96) 



and do perturbation on R{b) instead of doing perturbation directly on Dp^ib). This way the 
Borel transform can be better described in the most important region in the Borel integral, 
i.e., between the origin and the first IR renormalon singularity and the region just beyond 
the singularity ||21 . 



Nonetheless, it is interesting to observe that the estimated value (|95[) is remarkably close 
to the recent estimate of gluon condensate 0.006 ± 0.012 GeV^ which was obtained by 
fitting the spectral function of hadronic r— decay using QCD sum rule. 

Finally, some comments are in order. A full amplitude in QCD in general has infinitely 



many cut singularities in the complex coupling plane as shown by 't Hooft [^. One may 



wonder how this can be compatible with our proposed mechanism that is based on the 
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proposition that the nonperturbative amphtude as well as the Borel resummed amplitude 
of perturbation theory have cuts only along the positive real axes. The resolution of this 
question lies probably with the nonconvergence of the OPE. Since the OPE is expected 
to be an asymptotic expansion |3^, each term in the OPE, to the Wilson coefficient of 



which the Borel resummation is applied, does not have to have the same singularities of 
the true amplitude. Also, throughout this paper, the nonperturbative effects due to the 
higher dimensional operator condensates were consistently ignored. We expect, however, 
there should be no fundamental difficulty in incorporating them along the lines described in 
this section. We conjecture that the nonperturbative effect by the condensate of a dimension 
2n operator in the Adler function can be written as 

-DNp(as) = Cni-asY"'" exp(-n//3oas)C'2n(as) , (97) 

where f„ is a constant calculable from the RG equations on the Wilson coefficient and the 
condensate, and C2„ is a modified Wilson coefficient defined in a similar fashion as the C4 
in Eq. ([7^). Here C2n should be Borel resummed in the manner described in this section. 
The unknown constant C„ can then be determined by demanding that the imaginary part 
of order a^'^" exp(— n//?oQ!s) be canceled by the imaginary parts coming from the amplitudes 
associated with the operators of lower dimensions. 



VI. SUMMARY 

Based on the general argument on Borel resummation of a same sign perturbation series 
we have argued that the nonperturbative effect associated with the divergence of the per- 
turbation series should have a branch cut only along the positive real axis in the coupling 
plane. Demanding that the nonperturbative amplitude have such a branch cut constrains 
the form of a part of the nonperturbative amplitude, the part from which the branch cut 
arises, sufficiently that a relation can be established between the usually incalculable real 
part of the nonperturbative amplitude and the imaginary part that is calculable from the 
Borel resummation. This way part of the real part of the nonperturbative amplitude, which 
usually includes the leading term in weak coupling expansion, can be calculated from Borel 
resummation of the perturbation theory. 

As a nontrivial test, this mechanism was applied to the ground state energy of the dou- 
ble well potential and the two-point function in the two-dimensional 0{N) nonlinear sigma 
model at order In agreement with our proposed mechanism the nonperturbative am- 

plitudes in these models have branch cuts only along the positive real axis in the coupling 
plane. With this mechanism the leading terms of the nonperturbative amplitudes in these 
models could be calculated with good accuracy from the first terms of the corresponding 
perturbation series. 

We then applied this mechanism to the QCD condensate effects, particularly the gluon 
condensate effect, and suggested that some of the condensate effects can be calculated from 
perturbation theory, and gave an estimate of the nonperturbative amplitude induced by 
the gluon condensate using the known perturbative calculations of the Adler function. We 
observed that this mechanism could be applied to QCD, despite the fact that a true QCD 
amplitude has an infinite number of cut singularities in the coupling plane, since the OPE 
to which the Borel resummation is applied is not a convergent expansion. 
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